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Abstract 


In this paper, a bottom-up energy system optimisation model is proposed in order to support 
planning policies for promoting the use of renewable energy sources. A linear programming 
optimisation methodology based on the energy flow optimisation model (EFOM) is adopted, 
detailing the primary energy sources exploitation (including biomass, solid waste, process by- 
products), power and heat generation, emissions and end-use sectors. The modelling frame- 
work is enhanced in order to adapt the model to the characteristics and requirements of the 
region under investigation. In particular, a detailed description of the industrial cogeneration 
system, that turns out to be the more efficient and increasingly spread, is incorporated in the 
regional model. The optimisation process, alming to reduce environmental impact and econ- 
omical efforts, provides feasible generation settlements that take into account the installation 
of combined cycle power plants, wind power, solid-waste and biomass exploitation together 
with industrial combined heat and power (CHP) systems. The proposed methodology is applied 
to case of the Apulia region in the Southern Italy. 
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1. Introduction 


As widely acknowledged, energy consumption is one of the most reliable indi- 
cators of the development and quality of life reached by a country and the necessity 
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of satisfying a forecasted energy demand, over a certain time period, is the basis of 
energy planning. In particular, energy planning builds and verifies strategies in 
energy economy, which is, using the definition of World Energy Council, “that part 
of economics applied to energy problems, taking into account the analysis of energy 
supply and demand, as well as implementation of the means for ensuring coverage 
of energy needs in a national or international context” [1]. 

The energy planning discipline must take into account political aspects, social and 
environmental considerations, and is carried out taking into account the historical 
data collected in the previous energy plans of the country under examination. 

Energy planning methods are generally classified in three categories [2]: planning 
by models, by analogy and by inquiry. The accuracy of these methods depends on 
the time interval under investigation: short-term and medium-term (up to 10 and 20 
years), long-term (beyond 20 years). 

The planning by models methodology includes the econometric model and 
optimisation model. The econometric model generally relies on mathematical and 
statistical methods (such as regression analysis) to study economic systems. In parti- 
cular, its alm is the empirical testing (validation) of theoretical models, as well as 
the derivation of quantitative statements about the operation of economic aggregates 
[1,3]. All the econometric models are based on the use and implementation of statisti- 
cal data. They deal with problems implying one or several energy forms, different 
energy sectors and energy uses. Seldom do they take any explicit optimisation 
approach into consideration. Sometimes an implicit optimisation tendency is under- 
stood. 

The optimisation model, when the approach of the best possible solution according 
to a goal function is required, makes the step from a description by a model to a 
prescription by a model, as an optimisation procedure will demonstrate that any 
deviation from a determined situation leads to a degraded one [4,5]. This is the most 
important and broadest category of tools for energy planning. In particular, great 
relevance goes to the family of the multi-period linear programming models [6-8]. 
In the following, some of these models are described. The Brookhaven energy system 
optimisation model (BESOM ) attaches all costs to energy flows and minimises their 
sum over 1 year, while the time-stepped energy system optimisation model (TESOM) 
makes consecutive BESOM-type optimisations for single years [9]. The market allo- 
cation model (MARKAL), a successor of BESOM, is a large scale, technology- 
oriented activity analysis model which integrates the supply and end-use sectors of 
an economy, with emphasis on the description of energy related sub-sectors [10,11]. 
It can include 1 to 16 time periods of the same length, each one indexed with the 
central year. The multiple energy systems of Australia (MENSA) [12] is an improved 
and regionalized version of MARKAL that chooses the combination of demand-side 
and supply-side technologies which delivers the energy services at least cost, aver- 
aged over a specific time period. Analogously to MARKAL, the energy flow optimis- 
ation model (EFOM) provides an engineering oriented bottom-up model of a national 
energy system and has been developed under the approval of the Commission of the 
European Communities [13]. The EFOM describes the energy system as a network of 
energy flows, by combining the extraction of primary fuels, through a number of 
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conversion and transport technologies, to the demand for energy services or large 
energy consuming materials. In the EFOM, the planning horizon is defined by a 
certain number of periods, generally of different length, consisting of two seasons 
each divided into two periods [14,15]. A further application of the EFOM model, 
named energy planning optimisation model (EPOM), is adopted in [16], where the 
planning horizon is split into time periods of the same length each divided into no 
more than eight sub-periods. 

The planning by analogy methodology [17] allows the simulation of the same 
quantity, with a time lag, in a less developed country, through the use of a leading 
case as reference and the knowledge of the time behaviour of a quantity in a more 
developed country. The “analogue” approach is often used to check and compare 
outputs produced by other methods. 

The previously mentioned tools may loose reliability when a long-term investi- 
gation is carried out, although MARKAL model covered a 45-year span when 
involved in greenhouse gases (GHG) emission abatement optimisation [18]. In this 
case, the inquiry system, called Delphi method [19,20], is applied, based on question- 
naires submitted to a selected panel of experts, in order to achieve, by statistical 
evaluation of their answers, an accurate chart of the future. 

Independently of the method, the energy planning requires a detailed preliminary 
study of the energy system. This important step assumes that a careful observation 
of all the phenomena involved in the evolution of energy demand and supply has 
to be carried out. 

Recently, the energy systems are undergoing a development trend characterized 
by the following principal guidelines [21,22]: 


e The privatisation of the most important energy sectors (electricity and natural 
gas), turning the previous monopolies into free competition among different com- 
panies; in particular, the unbundling of vertically integrated energy companies 
occurred in the electricity sector, by splitting the generation, transmission and 
distribution activities; 

e The community growing awareness about the environmental impact caused by 
large conventional power plants, joined with a greater interest towards distributed 
generation technologies based upon renewable resources and cogeneration; 

e The energy planning activities as regional concerned instead of national. 


In such a scenario, the energy planner has to shift the border of the system under 
study towards a smaller observation area where, invariably, several new constraints 
of different nature are involved. 

This paper aims to provide a methodology for regional energy planning over a 
time interval of some decades. To this purpose, the typical modular structure of the 
EFOM, usually applied for supporting planning policy in a whole country [23,14], 
has been tailored to the requirements of a regional energy system, including the 
description of primary energy source exploitation, power and heat generation, emis- 
sions and end-use sectors. A detailed description of the industrial cogeneration sys- 
tem has been incorporated in the model so that a more realistic representation of the 
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actual industrial sector of the region under investigation is provided. The results from 
scenarios on strategies to reduce emissions and to improve the energy production 
from renevvables turn out to be useful for a real implementation on actual regional 
energy system vvith reasonable costs. 


2. Energy system structure and general assumptions 


The proposed optimisation method aims at the determination of an optimal mix 
of technologies for the energy system, subject to a number of boundary conditions, 
such as emission limits, cost reduction, etc. The overall structure of the adopted 
model is shown in Fig. 1. 

The energy system is represented as a network of energy chains, starting from the 
primary energy supply and ending in the end-use sectors. The model is driven by 
an exogenous demand for useful or final energy. Follovving the energy chain rep- 
resentation, the model is built, in a modular vvay, into sub-systems according to 
relevant region specific availabilities and requirements. The sectors incorporated in 
the model include primary supply sectors, the povver and heat generation sectors and 
end-user sectors: industrial, transport, agriculture and fishery, residential and com- 
mercial. 


2.1. Primary energy supply sector 


The primary supply sector includes fossil fuels (coal, oil, natural gas, etc.), indus- 
trial by-products (blast furnace gas and cokery gas) and local renewables (biomass, 
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Fig. 1. Modular structure of the regional energy system. 
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solid waste and natural resources). All the primary energy flows are evaluated in 
tons of oil equivalent (toe) units. This sector provides for electricity and heat pro- 
duction deriving from large-scale as vvell as small-scale technologies and for other 
needs directly in the end-use sector. 


2.2. Intermediate conversion of primary energy 


The large-scale technologies are adopted both for centralized electricity generation 
and for industrial povver and process steam production. VVhereas small-scale techno- 
logies largely dominate distributed electricity production by local renevvable exploi- 
tation and civil thermal/electric energy production. The vvhole net electricity pro- 
duction delivered to the povver grid comes from povver plants, distributed generation 
and from industrial combined heat and povver (CHP) plants. The total generation 
capacity to be installed takes into account the possibility to export energy to the 
neighbouring regions. 

Various technologies, such as conventional thermoelectric plants, combined eycle 
plants, photovoltaic panels, micro-turbine systems and etc., have been supposed for 
intermediate conversion of primary energy. The choice and the adoption of these 
technologies depends on the objective of the optimisation process: if a greater impor- 
tance is given to economical requirements, cheapest options vvill be considered: on 
the other hand, vvhen the goal is to reduce pollutant emissions, the environmentally 
friendly technologies are preferred. 

Due to the wide adoption of CHP facilities in industrial processes, a particular 
care is spent for the description of the relevant cogeneration system. On the basis 
of the sample cogeneration system schemes, reported in [24,25], the cumulative rep- 
resentation of the whole industrial cogeneration system producing high, medium and 
low pressure steam has been obtained, as shown in Fig. 2. 

The industrial cogeneration system includes gas turbines (GT's) and/or extraction 
condenser steam turbines (STs) for power generation. The STs may be fed by high 
pressure (HP) boilers and/or by the high pressure steam from heat recovery steam 
generators (HRSGs), which are set downstream the exhaust of GTs. The high press- 
ure boiler system may fire various fuels such as oil, blast furnace gas and cokery 
gas, while the combustion chamber (CC) of the GTs fires only natural gas. A cascade 
system of pressure reducing valves (PRVs) grants the availability of three process 
steam levels. 


2.3. End-use sector 


The end-use subsystem defines a set of energy demand disaggregated in electric 
power, industrial process steam, civil heat and other needs (coal for blast furnace, 
natural gas for cooking systems, petrol and diesel oil for means of transportation, 
etc.). In particular, civil and industrial users can satisfy their electricity needs both 
withdrawing energy by the grid and self-producing power. However, electric energy 
production by civil sector cannot be delivered to the grid. The other sectors are only 
allowed to buy electric energy from the grid. 
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Fig. 2. Schematic outline of the industrial cogeneration system. 


2.4. Environmental impact assessment 


Two main aspects concerning environmental impact are considered. The first 
aspect deals with the quantification of emissions to air, most notably CO», NO,, SO,, 
PM o, due to consumption of fuels. To this purpose, suitable specific emission fac- 
tors, able to assess the tons of pollutant per toe of burned fuel, are adopted. The 
pollutant emissions are limited below a certain tolerance threshold, such as the Kyoto 
Protocol (KP) requirements, by proper inequality constraints in the mathematical for- 
mulation. 

The further quantification of environmental impacts consists in evaluating the 
resulting damage costs following the ExternE methodology [26]. The ExternE pro- 
cedure is a bottom-up methodology assessing the whole life and fuel cycle of a 
specific plant. It employs an impact pathway approach tracing the emissions from 
the source to the impact, assessed by means of dose—response functions. A broad 
variety Of burdens is considered, starting from impacts on human life and health up 
to visual amenity. The monetary valuation of the resulting welfare losses follows 
the approach of “vvillingness-to-pay” for improved environmental quality [26]. In this 
paper, the whole regional “vvillingness-to-pay” is managed to fit with medium exter- 
nal costs and is included in the objective function to be minimised in the optimis- 
ation process. 
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2.5. Multi-period approach 


In order to simulate the dynamics of the energy system, over a period typical of 
the regional energy planning (10-20 years) the total study period is divided into sub- 
periods (years), and the solution of the model is obtained by optimising an objective 
function over the vvhole study period. The annual demand of electricity, industrial 
process steam and civil thermal power is represented by a load duration curve (LDC) 
that is fitted by a step-wise behaviour, as reported in Fig. 3. 

Each step corresponds to the demand level expected throughout the relevant time 
interval. Practical computational reasons allow a maximum of eight time steps for 
the load duration curve fitting. 


3. Model formulation 


The EFOM-based optimisation model, described in this section, aims to determine 
the optimal use of resources in energy supply and conversion and, to some extent, 
in final energy demand. The total present cost of the entire energy system is 
optimised by a linear programming procedure over the whole study period. The 
resulting objective function is minimised in the presence of suitable equality and 
inequality constraints. In this section, after a list of all used indices, parameters and 
variables, a detailed description of the cost function and of the equality/inequality 
constraints is given. 


3.1. Indices 


t represents mid-time periods (year): t = 1,....N;; 
p represents the time intervals of LDC: p = 1,..., N, and N, = 8; 
i stands for primary energy forms: 


demand 
level 


hours 8760 


Fig. 3. Step-wise curve fitting the annual LDC. 
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i — 1 fuel oil 


i — 2 coal 

i — 3 natural gas 
i — 4 diesel oil 

i — 5 petrol 


i — 6 liquid propane gas 

i — 7 municipal vvaste 

i — 8 biomass 

i — 9 blast furnace gas 

i — 10 cokery gas 

i = 11 = N, refinery gas; 

j represents the primary energy conversion options: 


centralized electric energy production plants 
j= 1 multi-fuelled thermoelectric 
j = 2 combined cycle 


distributed electric energy production plants 
j =3 hydroelectric 

j =4 photovoltaic 

j =5 wind 

j= 6 solid waste-to-energy 

j “ 7 biomass-to-energy 


civil thermal/electric energy production facilities 
j =8 solar thermal collectors 

j =9 photovoltaic roofs 

j= 10 gas micro-turbines 


industrial CHP facilities 
j= 11 gas turbines 
j= 12 = N, steam turbines; 


k stands for fuel mix inputs of thermoelectric plants: 
k= 1 oil (30%)—coal (70%) 

k = 2 oil (40%)—coal (60%) 

k = 3 = N, oil (50%)-coal (50%); 


s stands for end-use sectors: 

s = 1 residential and services 
s = 2 industry 

s = 3 agriculture and fishery 

s = 4 N, transports; 


r represents the atmospheric pollutants: 
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1 CO, 
2 NO, 
=3 SO, 
4=N, Particulates. 


3.2. Parameters 


Generally speaking, a value which remains constant in a study or for which only 
two or three specifically fixed amounts are assigned, is considered as a parameter. 
In our model, the following parameters are adopted: 


T,, Operation time fraction of the LDC p-th interval (h) 

6 corrective factor taking into account a reserve margin on annual load peak 
o: energy transmission line losses factor 

DR discount rate 

€;, r-th pollution factor per unit of the i-th fuel consumption (ton/toe) 


Energy conversion system 

AV, availability factor of the /-th installation type (h/y) 

A, amortization period of an investment in the /-th installation type (y) 
D, life period of the /-th installation type (y) 

Hi, mass percentage of the i-th fuel in the k-th mix 

c;; unit consumption of the i-th fuel by the j-th conversion option (toe/MWh) 
Industrial sub-system 

h, unit consumption of the i-th fuel by HP steam boilers (toe/ton) 

m, unit consumption of the i-th fuel by MP steam boilers (toe/ton) 

Ocr linear gas turbine power-to-exhaust flow curve slope (MWh/ton) 

Or, Op, Ow linear steam turbine power equation coefficients (MWh/ton) 

v HRSG specific exhaust steam outflow rate (ton/ton) 

Am Am AL PRV efficiencies 

y medium unit extraction limit (ton,,,/ton,,) for steam turbines 


Civil sub-system 

1), i-th primary energy conversion efficiency of boiler systems 

hpr, heat-to-power ratio (MW/MW.) of gas micro-turbines at the p-th load con- 
dition. 


3.3. Variables 


The choice of variables is based on the degree of freedom and flexibility that the 
planner means to attribute to some quantities. In this case, a distinction between 
“explanatory” and “explained” variables is adopted [2]. 
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3.3.1. Explanatory variables 
The “explanatory” variables generally define the evolutionary scenarios of the sys- 
tem. For our purpose, the following “explanatory” variables are used: 


Fi.” Fi. lower and upper bounds for the i-th fuel consumption at the t-th year (toe) 
RP,, available energy potential for renevvable-based generation option (j = 3-9) at 
the f-th year (toe) 
B,, total capacity of the j-th generation option, already existing at the beginning of 
the planning horizon and still in operation at the t-th year (MW) 
Exp, electric power exported during the p-th interval of the th year (MW) 
PP,,, purchased electric power by the s-th end-use sector during the p-th interval 
of the ¢-th year (MW) 
electric power amount delivered to the grid and self-produced by the s-th end- 
use sector (SP, ,,F0 only for s = 2), during the p-th interval of the ¢-th 
year (MW) 
OP.» electric power demand of the s-th end-user during the p-th interval of the t- 
th year (MW) 
OT, civil thermal power demand at the p-th interval of the t-th year (MW) 
®,, share of civil thermal demand in the 1-th year satisfied by firing the i-th energy 
form (%) 
QH,» QM,» OL, , high, medium and low pressure industrial process steam demand 
during the p-th interval of the ¢-th year (ton/h) 
QF,,, demands for the i-th fuel by the s-th sector in the f-th year (toe) not included 
in the optimisation process 
TOL,, emission limits of the r-th air pollutant in the f-th year (ton). 


SP 


sp 


ep 


The set of the “explanatory” variables also include the unit costs used to assess 
annuities in the objective function. The description of these further variables will be 
provided in Section 3.4. 


3.3.2. Explained variables 

In each planning study, one or more “explained” variables have to be evaluated 
within a logical framework. The level of detail, chosen for the proposed method- 
ology, suggests the adoption of the following 'explained' variables: 


T, additional capacity of the j-th installation type at the t-th year (MW) 

R;,, recovered thermal power from the j-th generation option during the p-th interval 
of the t-th year (MW) 

P;,,, power output of generation option j = 1, from the k-th mix, during the p-th 

interval of the ¢-th year (MW) 

P,,, power output of the j-th generation option during the p-th interval of the t-th 

year (MW) 

industrial HP steam boilers outlet from i-th fuel during the p-th interval of the 

t-th year (ton/h) 


H 


itp 
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Sip Steam turbine inlets from the HP boilers burning the i-th fuel (i = 1,9,10,11) 
during the p-th interval of the t-th year (ton/h) 
S,, steam turbine inlets from HRSGs during the p-th interval of the 1-th year (ton/h) 


Vip» Vip steam flows degraded by PRV during the p-th interval of the t-th year (ton/h) 
M;,, industrial MP steam boilers outlet from i-th fuel during the p-th interval of the 
t-th year (ton/h) 


E, ,E,, extractions from steam turbines during the p-th interval of the r-th year 
(ton/h) 
W,, condensing flow from steam turbines during the p-th interval of the ¢-th year 


(ton/h). 


The main part of the above defined “explained” variables refers to the industrial 
cogeneration model detailed in Appendix A. 


3.4. Objective function 
The objective function consists in the total actualized cost C7 of the primary energy 
conversion over the selected time horizon: 
Cr = CI + CF + CV + CE (1) 


where CI is the total actualized investment cost of the installed plants, CF represents 
the total fixed cost, CV takes into account the variable costs and CE the external costs. 
The investment costs CI is composed of two terms: 


CI = CIC + CP (2) 


CIC being the capital investment cost and CP? the stripping cost. Each of these terms 
is obtained as sum of annuities over the period of plant lifetime and assumes the 
following expression: 


= 2 a + uam 5 KI, xl; (3) 


2” 


where KT, : is the unit capital cost (€/MVV) of the j-th installation type at the 7-th 
year, Dam, = {max (Lt - Aj + 1),...,t } and AF; the j-th investment type amortiz- 
ation factor: 


p-o 
an= Sufi Sy 


p being the interest rate of capital and 6 the inflaction rate. 
The stripping cost CP? is expressed as follows: 


cr = 5 SE a + AU (4) 


JeOsrrpt = 1 
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where KS; n; is the unit stripping cost (€/MVV) of the j-th generation option installed 
at the t-th year, Qsrrp is the set of all je (1,...Nj] ə” N, - Dj > 0. 
The total actualized fixed cost CF has the follovving expression: 
N; 1 Nj 
t=1 j=l teOpuR,, 


where KFX,, is the unit operation and maintenance cost (€/MW) of the j-th instal- 
lation type at the 7-th year and (lpuz,, = {max(1,-D,; + 1).....t]. 

The cost CV represents the procurement costs for all fuels necessary for the electri- 
cal energy generation and process steam production over the planning horizon: 


N, 
. 1 
CV — 2 a+ pay aug] Zam (Səs X Pu? m 


kel 1 p=1 


+ $ XP Py 1,) + > KFL, ch; S Het T, (6) 
j= 


j=2 p=1 ieOnHpB p=1 
Np 
+ >, KFL,-m, ` Mp Tp) 
ieOmpB p=1 


where KFL,, is the unit price (€/toe) of the i-th fuel in the th year and Oypg = 
Omps = {1,9,10,11}. The coefficient U; representing the percentage of the i-th fuel 
in the thermal multi-fuel units, weights the produced power [27]. 

The external costs consist in the economical estimation of burdens occurred to 
people and environment because of energy chains including the life-cycle both of 
primary energy sources CEFand power generation plants CEP: 


CE = CE" + CE" (7) 
For the cost CE” the following expression is adopted: 


N 


CEr = 2 a ETT 23 er, br -T A 
* 5, cu Ba 1,) + x KEP, A YS Hus Tp (8) 
f-2 p=1 icÜmpB p=1 
Np 
+ M KEF,m, b Mis ni 
ie OmPpB p=1 


where KEF), is the unit cost (€/toe) assessing impacts related to transportation, treat- 
ment and utilization of the i-th fuel in the t-th year. 
As well as for the cost CE” the following expression is assumed: 
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>- N, 
x KEL l; vi KES; nn; 
M 2 2 (L * DRy ^ 2 be ( 7: e 
j=1t= jeQsTrpt= 1 
N, 1 N; 


TEQ 
S DUR; 


where KEI „ KEL, ,, KES,, (€/MVV) take into account the impacts related respectively 
to installation, operative life and stripping of the j-th generation option in the t- 


th year. 
3.5. Equality and inequality constraints 


The optimisation procedure has to run within a logical framework that drives the 
evaluation of the “explained” variables. To this purpose, the definition of the follow- 
ing relationships is needed. 


3.5.1. Construction time constraints 

At the beginning of the time horizon, construction time is taken into account, by 
adding in the optimisation procedure suitable constraints able to delay the energy 
production of the new installations. The time delay depends on the particular type 
of plant/facility installation. As a consequence, the following equality constraints 
are included: 


labs 1—07—1 
La = Lo = J = 2,6,7 (10) 
Li = 0 j = 3,4,5,11,12 


For the generation options corresponding to j = 8,9,10 it is assumed a construction 
time less than 12 months; thus the installation of this options can occur starting from 
t=. 


3.5.2. Peak demand satisfaction 

The first time step of each defined LDC corresponds to the forecasted annual peak 
load. Then the total electricity generation capacity in this period has to cover, with 
a suitable reserve margin, the internal load demand and the exported power: 


Nj Ny 

> İb. + > MET + əl Y Pu: + Esp) t (11) 
j=l TeOpuR,, s] 
j*38 
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3.5.3. Plant/facility operation limits 
The net power output of each generation option cannot exceed the relevant 
installed capacity: 


»n Pico B; + > La j= 


k=1 TeOpun,, (12) 
P, SB, + py Ii. J = 2.--..12 
TEODUR, 


for 1,...,N, and p=1,....N,,. 


3.5.4. Limits on renewable energy potentials 

Limits on energy production (heat and/or power) from renewable and local energy 
sources have to be considered. Limitations are imposed directly on the power output 
of the generation options: hydroelectric, photovoltaic power plants, wind, solar ther- 
mal and dam roofs. The corresponding constraints are expressed as follows: 


T iss y Pi; T,ERP,, j = 3,4,5,8,9 (13) 
p=1 
where 4.86 is the gross MWh-to-toe ratio, for t = 1,...,N, and p = 1,..,N,. Limits 
on energy production by local renewable fuels, such as biomass and municipal waste, 
are taken into account by fixing a maximum consumption threshold and are illus- 
trated in Section 3.5.10. 


3.5.5. Limits on electrical energy generation 

Through the availability factor AV, (hours/year), maintenance and failure periods 
are defined and, contemporaneously, the maximum energy production is obtained 
for each generation option. Therefore, the annual energy production does not have 
to exceed this threshold, that is: 


No 7 Ne 
> ( 2 Phase) TSA (By de > » j=l 
p=1k=1 TEQDUR,, 
i (14) 
Y Pir 1,5 AV; İb. + > b) j = 2,...,12 
p=] TEQDUR,, 
for t = 1,...,N,. 


3.5.6. Electricity generation and consumption balancing 
The balance among the electricity generation, the internal consumption and the 
net povver exchange vvith the grid has to be satisfied. The follovving expressions 
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establish the electric povver balance for each economic sector and for the overall 
regional system during the planning horizon: 


2555” (15) 
j= 940 
2000002851 00: (16) 
j=11,12 

PP. > OPs 5 (17) 
PPR az ES QP,,, (18) 
Ny 7 N, 

y Pi + Pup + SPo tp E ce exp, + py PP] (19) 
k=1 j=2 s=1 


for t = 1,...,N, and p = 1,...,N,. It should be noted that the produced power in the 
civil sector is assumed to cover the sole internal needs. Whereas the industrial sector 
is allowed to deliver energy to the network. 


3.5.7. Industrial CHP facilities 

The most common cogeneration options used for industrial applications are: simple 
and double-extraction steam turbines, back pressure steam turbines, gas turbines with 
waste heat boiler and combined cycles [28]. As shown in Fig. 2, the purpose of our 
model is to include all these options in a simple scheme. 


3.5.7.1. Inlet boundaries The steam turbine system may receive input steam both 
from HP boilers and HRSG systems. These input flows derive from leading flows, 
namely the whole output of HP boilers and HRSG systems, that cannot be exceeded: 


Sip Hp Vie Oypg (20) 
V : 
== Pu j=11 Q1) 
OGT 


for t = 1,...,N, and p = 1,...,N,. 


3.5.7.2. Power balance Assuming linear steam-flow/electric-power characteristics 
for the steam turbine generators and negligible heat loss, the power equation may 
be written in the form: 


Pity = Of E,, + OgE,, t Op W,, j= 12 (22) 


Lp 


for t = 1,....N, and p = 1,....N,. 


3.5.7.3. Mass balance constraint The sum of the input steam flows of the turbine 
system has to be equal to the sum of the output steam flows: 
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Sip + > Sirp - E + E: + ye (23) 


¡SOHPB 


for t = 1,...,N, and p = 1,...,N,. 


3.5.7.4. Extraction boundary The sum of the extraction flows is constrained below 
a fixed portion of the total input: 


Esp + B,erİs, ag BS So) (24) 


ieOkhpB 


for t = 1,...,N, and p = 1,...,N,,. 


3.5.6. Industrial process-steam balance 

Since the self-produced electric power is already fixed in Eq. (16), the thermal 
load following is carried out by integrating the HRSGs with boiler outputs. Referring 
to Fig. Al of Appendix A, the industrial process-steam balance is expressed by the 
following equations: 


P x. 0H.» + Vip 


1 . Jp = um 
2p Asp Suy) + VG Sew A şen (25) 
¡sQupB 
= Mi. + E,» + Aw Vip OM.» + Vip (26) 
ie QMPB 
Esp T Ar Vip = OL. p (27) 


for t = 1,...,N, and p = 1,...,N,. 


3.5.9. Civil thermal power non-conventional production 

In order to satisfy hot water and heating requirements, civil users may either use 
conventional boilers or non-conventional facilities as solar collectors and micro-tur- 
bines. The thermal power, recovered from the non-conventional facilities, has to 
respect the following conditions: 


Ris» — Pi J=8 (28) 
Rip hpr Pi, j=10 (29) 
Y. RP. (30) 


j = 8,10 


for t = 1,....N, and p = 1,...,N,. Eq. (30) imposes that the total recovered power 
cannot exceed the internal needs. 


3.5.10. Primary energy consumption 
Fuel consumption is related to plant/facilities power production, conventional civil 
thermal plants output, industrial boiler output and other needs. The whole sum of 
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these contributions, for the i-th fuel at the /-th year, has to be comprised between 
suitable lovver and upper boundaries: 


11 
Lc Y cta Pu T, + 351 T 


k=1 p= j=2 p=1 
Np 


Np 
— 12.44. var. (or. 2. MES + +h b, T, (31) 


j = 8,10 p=1 


T Mi; x I, onun, 

p=1 
where h; = O(m, = TO and 12.44 is the net MWh-to-toe ratio; for 
—],...,N; and t=1,....N,. 


3.5.11. Environmental constraints 

The unit consumption of an energy form is related to pollutants emission through 
suitable emission factors. Therefore, the total annual emission of the r-th pollutant 
is evaluated with the same approach adopted for fuel consumption assessment. The 
environmental scenario defines the time evolution of the emission tolerance to be 
satisfied by the annual emission level of each pollutant: 


N, Ny 
ea T ` €, Y PaT, + 


i=l k=1 p=1 j=2 p=1 


P 


N; 
ai 1, 2 (or..- 2 Bəkr, +2 OF (32) 
= j= 8,10 =i 
Np Np 
+ Y erh) HaT, >, ewm 2, Mu, T,STOL,, 
ieOnpB p=1 ieOmpB p=1 


for r = 1,...,N, and t = 1,...,N,. 


4. The energy system under study 


The procedure described is applied to the assessment of optimal energy plans for 
Apulia region, in southern Italy. The actual snapshot of the energy system is assumed 
as the starting point of the planning horizon. Therefore, following the modular struc- 
ture illustrated in Fig. 1, the primary energy supply sector, the intermediate energy 
sector, the end-use demand (electricity and heat) sector and the environmental impact 
are detailed. 

Due to the presence of great iron and steel production sites, a remarkable quantity 
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Table 1 
Energy balance of the system under study 


Solid Oil Gas fuels Renevvable Electric Total 
fuels products (ktoe) energy (ktoe) 
(ktoe) (ktoes) (ktoe) 
Fuels Non- 
(ktoe) fuels 
(ktoe) 
Gross consumption 6528 4940 1570 111 89 -958 12 280 
Electric energy -3697 -1252 0 -111 -89 5149 0 
conversion 
Consumptions/losses of -893 152 -20 0 0 -2797 -3558 
energy sector 
Other uses 128 815 0 0 0 943 
End-use consumption 1810 3025 1550 0 0 1394 7779 


of by-product gases is available and utilized in industrial energy self-production. 
Moreover, Apulia imports from neighbouring regions and countries the greatest part 
of the fossil fuels. 

The forecasted gross energy balance of the system under investigation is reported 
in Table 1. The data, illustrated in this table, come from a previous study developed 
in Ref. 1291. The primary energy sources are classified as solid fuels (coal, blast 
furnace gas and cokery gas); oil products (fuel oil, diesel oil, petrol, LPG and refinery 
gas); gas fuels (natural gas); renevvable fuels (municipal waste, biomass) and renew- 
able non-fuels (hydro, wind, photovoltaic). 

It is worth remarking that the actually exploited renewable fuels are biomass 
whereas the main electric energy from renewable non-fuels is produced by wind 
power plants. 

In Table 2, the installed electric generation capacity is reported for every existing 
generation option. The total generation capacity of the region is 4850 MW. 

The main part of the conventional thermoelectric capacity (2447 MW) has been 
operating since 1998 whereas the remaining 625 MW since 1980. The hydro-electric 
plants started working in 1987, as well as the photovoltaic plants. The wind power 


Table 2 
The installed electric generation capacity 


Generation option Installed capacity at 2002 (MW) 
Thermal power plant 3102 

Hydro power 2 

PV power 0.6 

Wind power 163 

Biomass power 43 


Industrial steam turbine 1540 
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Table 3 
The End-use consumption 


Civil Industry Agr and Transport Total(ktoe) 
fishery 
Solid fuels (ktoe 937 1776 0 0 1809.7 
Oil-Products (ktoe) 181.6 470.6 328.9 2044.3 3025.4 
Gas-Fuels (ktoe) 692.2 853.6 2.4 2.0 1550.2 
Electric Energy (ktoe) 584.9 745.6 50.1 13.3 1393.9 
Total (ktoe) 1492.5 3846.0 381.4 2059.6 7779.0 
Table 4 
Emission to air total amounts 
CO, (ton) NO, (ton) SO, (ton) Particulates (ton) 
Solid fuels 7 976 457 23 591 35 145 25 898 
Oil products 8 971 610 78 972 22 976 6647 
Gas fuels 3 627 433 6356 0 310 
Electric energy 21 016 220 63 149 118 140 60 255 
Total 41 591 720 172 068 176 261 93 110 


plants have been installed in the period 1996-1998, while the biomass exploitation 
dates back to 1997. The industrial steam turbine installations were undertaken from 


1990-1997, though most of the installations (984 MW) occurred in 1997. 


The end-use consumptions of the region under investigation are illustrated in Table 


3, for every economic sector. 


It can be noted that thanks to the widely spread gas network, both the civil and 
industrial sectors greatly exploit natural gas for all purposes, whereas the oil products 
are essentially used for transport. 

The total emission amount of CO», NO,, SO, and particulates due to the end-use 
consumption of fuels, as reported in Tables 1 and 3, and electric energy production 


at year zero are reported in Table 4. 


Table 5 


Annual growing rates of the electricity and heat demand 


Civil Industry Agr and fishery Transports 
Electricity demand 1.10 1.93 1.95 0.30 
annual rate (76) 
Heat demand annual rate. 1.75 1.00 = = 


(%) 
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Fig. 4. Scenario 1: Generation capacity to be installed. 


5. Numerical results 


The proposed energy planning methodology is applied to the regional system illus- 
trated in Section 4. To this purpose suitable scenarios are considered taking into 
account different regional economic and environmental policies. 
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Fig.5. Scenario 1: Civil primary energy gross consumption for heating purposes. 
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Fig. 6. Scenario 1: Industrial primary energy gross consumption for steam production. 


5.1. Scenario 1 


The annual rates of the electricity required by every end-use sector, as well as 
the heat demand, are supposed growing as described in Table 5. 

It can be noted that, since no heat power is demanded by the sectors other than 
civil and industrial, no per cent increase is assumed for these sectors. The allowable 
emission amount of CO, is considered in accordance with the Kyoto Protocol 
requirements, i.e. 5-7% decreasing total rate, with regard to 1990 emission amount, 
for the first 10 years and then constant. Analogous limitations are assumed for the 
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Fig. 7. Scenario 1: Primary energy gross consumption for electricity production. 


other pollutants. Furthermore, all the unit costs, defined in Section 3, are supposed 
increasing with a rate able to offset the inflation rate. The external costs of the 
electricity production are neglected in this scenario. 

The resulting capacity to be installed over the planning horizon is reported in Fig. 
4 and Table 6. 


The main new production capacity, for the first 6 years and for the last 10 years, 


KP constraint 


CO, Emissions [Mton] 
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Fig. 8. Scenario 1: CO, emissions to air. 
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Fig. 9. Scenario 1: Other pollutant emissions to air. 


derives from combined cycle power plants which turn out to be the most economic 
and efficient technology, whereas a great number of industrial steam turbines are 
installed at the 15th year because of the end-life of the industrial thermoelectric 
power plants. 

The primary energy gross consumption concerning civil heating is shown in Fig. 
5. It can be noted a growing use of the gas-fuels and a light decreasing use of the 
oil-products over the studied interval. 

The industrial process steam production requires the sole consumption of solid 
fuels and oil products as shown in Fig. 6. A reduction in the oil product consumption 
is noted in the first years in an attempt to minimise the costs, though the increased 
use of less expensive solid fuels (such blast furnace and cokery gas) yields higher 
emissions. This explains the turnround in consumption behaviour in the last 12 years. 

The gross consumption of the primary energy sources, due to electric energy pro- 
duction is reported in Fig. 7. In accordance with the policy of promoting efficient 
technologies and renewables, the growing installation of combined cycle power 
plants is pointed out by the gas-fuel gross consumption rise and the fall in solid fuel 
gross consumption. 

It can be further remarked that the electric energy from the renewable fuels and 
non-fuels experiences a small increase. 

The corresponding emission to air levels, over the time horizon, respect the fixed 
limitations. In particular, the CO, emission strictly follows the Kyoto Protocol con- 
straints, as reported in Fig. 8. 

In Fig. 9, the behaviour of the NO,, SO, and Particulates emissions is shown too. 
It has been verified that these emission values are well above the imposed limits. 
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Fig. 10. Scenario 2. Generation capacity to be installed. 


5.2. Scenario 2 


Assuming the same conditions of the previous case, the nevv scenario includes the 
external costs in the objective function. The external costs, related to MW of installed 
plant/facility and to toe of burned fuel, are obtained by adapting the ExternE method- 
ology study cases to the Apulia region situation [30]. The further limitation on 
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Fig. 11. Scenario 2: Industrial primary energy gross consumption for steam production. 


Time [years] 


Fig. 12. Scenario 2: Primary energy gross consumption for electricity production. 


environmental impacts gives rise to a new settlement of the generation capacity to 
be installed over the planning period as reported in Fig. 10 and detailed in Table 7. 

The main part of the combined cycle power plants (1462 MW) is installed at the 
3rd year and a remarkable installation (240 MW) of wind power can be observed 
at 11th year, whereas the installation of the industrial steam turbines at 15th year is 
still the same. 

The primary energy gross consumption for civil heating does not change, whereas 
a different behaviour can be observed for the solid fuel and oil product gross con- 
sumption for industrial steam production as shown in Fig. 11. 
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Fig. 13. Scenario 2: CO, emissions to air. 
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Fig. 14. Scenario 2: Other pollutant emissions to air. 


—4— scenario #1 —8— scenario #2 
7000 


6500 


6000 


5500 


5000 


4500 


4000 


0 12 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 


time [years] 


Fig. 15. Comparison between the total net capacities to be installed. 
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After a considerable decreasing in the first 2 years, the burned oil product assumes 
an increasing trend. On the other hand, the solid fuel consumption decreases in the 
first 3 years and then it keeps an almost constant value. 

The gross consumption of the various fuels, due to the electricity production, 
remarkably changes in the first years of the planning horizon because of the further 
restrictions on the environmental impacts, as reported in Fig. 12. 

The related emissions to air change as vvell, so that the strict dependence of the 
CO, and SO, emission on the solid fuel consumption is put in evidence. In Fig. 13, 
the CO, emission behaviour is illustrated and the other pollutant emissions are 
reported in Fig. 14. 

It can be noted that the inclusion of the external costs in the objective function 
drives the CO, emissions well below the Kyoto Protocol constraints. 

In order to summarize the main differences of the two scenarios, the total net 
generation capacity is drawn in Fig. 15. For each year of the time period, the net 
capacity includes the plants still in service and the new installations selected by the 
optimisation process. 

It can be noted an increase (about 800 MW) of the net generation capacity at the 
3rd year that is ascribable to the rise of the combined cycle installations observed 
in Fig. 10. A further deviation (about 200 MW) of scenario 2 from the scenario 1 
occurs starting form the 11th year of the planning horizon, owing to the greater 
employment of generation options such as wind power plants, waste-to-energy sys- 
tems, gas turbines for industrial cogeneration systems. 


6. Conclusions 


A bottom-up energy system optimisation model has been adopted to support 
regional energy planning policies. The proposed approach is based on EFOM modu- 
lar structure, detailing the primary energy sources exploitation (including biomass, 
solid waste, process by-product), power and heat generation, emissions and end-use 
sectors. Moreover, particular care has been given to the description of the industrial 
cogeneration system scheme that has been incorporated in the regional energy model. 
The effectiveness of the proposed methodology has been shown by carrying out 
simulations on the Apulia region energy system. Test results proved that the regional 
policy, almed satisfying the increasing heat and power demand by various end-use 
sectors through environmental friendly technologies, can be supported mainly by 
combined cycle installations and with less effort of wind to power, waste-to-energy 
and biomass exploitation and industrial cogeneration systems. The inclusion of the 
external costs, in the objective function permitted to further force the conventional 
thermoelectric power plants out of the energy planning, increasing the adoption of 
more efficient and less polluting facilities, such as cogeneration systems and wind 
power. Although, it needs to be remarked that other renewable options such as PV 
power turn out to be still not competitive because of the high capital costs. 
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Appendix A 


This appendix provides a detailed description of the industrial energy system, including 
cogeneration facilities already shown in Fig. 2 and all the external and internal linkages. 
The CC of the GTs fires only natural gas (i = 3), the GT outputs are both the electric 
povver, delivered to the utility grid or directly satisfying the internal demand, and the 
recovered HP-steam vvhich feeds the STs and satisfies the HP process-steam require- 
ments. The HP and MP boiler systems burn oil (i — 1), and all industrial by-products 
(i = 9, 10, 11). Their outputs supply steam to industrial processes at various pressure 
levels through a system of PRVs. The STs may be fed by high pressure boiler systems 
and/or the high pressure steam from HRSG systems. The ST outputs are the electric 
power delivered to the utility grid or withdrawn by industrial sector and the extracted 
steam suitable for satisfying the MP/LP process steam requirements. 
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Fig. A1 Overall outline of the industrial sector. 
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